Families of particles with different masses in PT-symmetric quantum field theory 



Carl M. Benderf] and S. P. Klevanskyff] 

"Physics Department, Washington University, St. Louis, MO 63130, USA and 

b Institut fur Theoretische Physik, Universitdt Heidelberg, Philosophenweg 19, 69120 Heidelberg, Germany 

(Dated: July 6, 2010) 

An elementary field-theoretic mechanism is proposed that allows one Lagrangian to describe a 
family of particles having different masses but otherwise similar physical properties. The mechanism 
relies on the observation that the Dyson- Schwinger equations derived from a Lagrangian can have 
many different but equally valid solutions. Nonunique solutions to the Dyson- Schwinger equations 
arise when the functional integral for the Green's functions of the quantum field theory converges 
in different pairs of Stokes' wedges in complex field space, and the solutions are physically viable if 
the pairs of Stokes' wedges are VT symmetric. 
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The standard model of elementary particles has three 
generations of fermions (leptons and quarks) whose 
masses range over several orders of magnitude. It is not 
known why there are three generations of masses and 
whether there are only three. This paper proposes a field- 
theoretic mechanism that might explain the occurrence of 
generations of particles having different masses but other- 
wise similar physical properties: There might be just one 
Lagrangian (or Hamiltonian) to account for the proper- 
ties of all these particles, but the functional integral con- 
structed from this Lagrangian may have many different 
physical realizations depending on the boundary condi- 
tions on the path of integration in complex field space. 
While the Dyson-Schwinger equations constructed from 
the functional integral are unique, the solution to these 
equations is not unique. The number of distinct solu- 
tions to the Dyson-Schwinger equations equals the num- 
ber of pairs of complex Stokes' wedges in function space 
in which the boundary conditions on the functional inte- 
gration can be imposed. For each pair of Stokes' wedges 
there corresponds a different field theory. 

Z. Guralnik et al pQ first recognized that for func- 
tional integrals, inequivalent classes of contours associ- 
ated with different complex boundary conditions give 
rise to nonunique solutions to the Dyson-Schwinger equa- 
tions. They argued that multiple solutions might account 
for inequivalent vacua. The key point of the current 
paper is that the pairs of Stokes' wedges in which the in- 
tegration contours terminate must be oriented in a PT- 
symmetric fashion in complex field space. If this is the 
case, there is strong evidence that the corresponding field 
theory will be physically acceptable; that is, the masses 
(poles of the Green's functions) will be real and the the- 
ory will be unitary. The mechanism proposed here is 
field-theoretic, but its application is not restricted to ele- 
mentary particle physics. Experiments on 'PT-symmetric 
optical wave guides [21 [3] and on PT-symmetric diffusion 
[3] have been reported recently. 



The conjecture discussed in this paper stems from re- 
cent research on VT quantum mechanics, where it has 
been shown that the 'PT-symmetric Hamiltonians 

H = p 2 + q 2 (iqy (e>0) (1) 

all have real positive spectra [SJ [5] . Each of these Hamil- 
tonians defines a conventional quantum theory with a 
Hilbert space having a positive inner product |7]. The 
time-evolution operator U — e~ lHt is unitary and thus 
probability is conserved. Spectral reality and unitary 
time evolution are essential for any quantum theory. 
These features are guaranteed if H is Dirac Hermitian. 
(By Dirac Hermitian we mean that H = W , where f rep- 
resents combined complex conjugation and matrix trans- 
position.) However, it is not necessary for H to be Dirac 
Hermitian for the spectrum to be real and for time evo- 
lution to be unitary; non-Dirac-Hermitian Hamiltonians 
can also define physically acceptable quantum theories. 

The Hamiltonians are VT symmetric because they 
are invariant under combined spatial reflection V and 
time reversal T ■ Such Hamiltonians are physically ac- 
ceptable because they are selfadjoint, not with respect 
to the Dirac adjoint f , but rather with respect to CVT 
conjugation, where C is a linear operator that represents 
a hidden reflection symmetry of H . The CVT adjoint 
defines a positive-definite Hilbert space norm. Not ev- 
ery PT-symmetric Hamiltonian has an entirely real spec- 
trum, but the spectrum is entirely real if and only if a lin- 
ear PT-symmetric operator C exists that obeys three si- 
multaneous algebraic equations [7]: C 2 = 1, [C,VT] = 0, 
[C,H] = 0. When the C operator exists, we say that the 
VT symmetry of H is unbroken. Finding the C oper- 
ator is the crucial step in showing that time evolution 
for a non-Hermitian PT-symmetric Hamiltonian is uni- 
tary. The phase transition between broken and unbroken 
regions for some PT-symmetric Hamiltonians has been 
observed experimentally 0I3]- 

The Hamiltonians in (fll) are smooth extensions in the 



2 



parameter e of the Dirac-Hermitian harmonic oscillator 
Hamiltonian (at e = 0) into the complex non-Hermitian 
domain (e > 0). As e increases from 0, the Stokes' wedges 
in the complex- 2 plane inside of which the boundary con- 
ditions for the eigenvalue problem 



ip"(x) + x 2 {ixY^(x) = Etp(x) 



(2) 



are imposed, rotate downward and become thinner. As 
shown in Ref. [5], at e = the Stokes' wedges are centered 
about the positive- and negative-real axes and have angu- 
lar opening 90°. At e = 2 the Stokes' wedges are adjacent 
to and below the real axes and have angular opening 60° . 
When e > 2, these wedges lie below the real axis. 

To illustrate the idea of this paper in a quantum- 
mechanical context we set e = 4 in (JlJ). The resulting 
x 6 Hamiltonian describes two different quantum theories 
because the eigenfunctions ip{x) can satisfy two differ- 
ent sets of boundary conditions [8] : (i) the conventional 
Dirac-Hermitian quantum theory for which ip{x) vanishes 
as | a; | — > oo in the complex- x plane in 45° wedges cen- 
tered about the real axes; or (ii) the unconventional VT 
theory, which is the extension in e of the harmonic oscil- 
lator. For this non-Hermitian quantum theory ip(x) also 
vanishes as |a;| — > oo in the complex plane in 45° wedges, 
but now these wedges are centered about argx = —45° 
and arga; = —135°. The one-point Green's function 
Gi = (x) distinguishes between these two theories. The 
conventional Dirac-Hermitian theory has parity symme- 
try, and thus G\ vanishes. The boundary conditions for 
the VT quantum theory violate parity symmetry, and as 
a result G\ has a negative-imaginary value. The nonva- 
nishing of G\ in the VI" theory is a purely nonpertur- 
bative effect; one cannot express G\ for the Hamiltonian 
H = p 2 + x 2 + gx 6 as a series in powers of g. 

The idea that different boundary conditions allow one 
Hamiltonian (or Lagrangian) to describe several different 
physical theories is general and extends beyond quantum 
mechanics to quantum field theories of fermion and/or 
boson fields of any spin and in any space-time dimen- 
sion. However, for brevity we consider here the massless 
Z?-dimensional pseudoscalar field theory (D < 2) hav- 
ing a selfinteraction of the form <j) in+2 (n = 1, 2, 3, . . .). 
(Under parity reflection <f> — > —<fi.) The Euclidean func- 
tional integral for the vacuum persistence functional in 
the presence of an external source J is 



Z[J] = (0|0) = f V(j)e 
Jc 

s = 



9 j4n+2 



4n + 2 



J(j> 



(3) 



At D = 1 , this quantum field theory reduces to a quasi- 
exactly-solvable quantum-mechanical theory [9]. 

For each integer n there are n + 1 different physical 
realizations of the quantum field theory in ([3| . To explain 
this we consider the analogous one-dimensional integral 



J c dip exp (— Lp 4n+2 ^ . When n — this integral exists 
only if the integration contour C begins and ends in the 
Stokes' wedges of angular opening 90° centered about the 
real-v? axis. These Stokes' wedges are shown in Fig. [T] 
The contour C must begin and end in different Stokes' 
wedges; if C begins and ends in the same Stokes' wedge, 
the integral vanishes. When n = 1, there are two possible 
choices for integration contour C; C may connect the two 
30°-Stokes' wedges centered about the real axis or C may 
connect the 30°-Stokes' wedge centered about —120° to 
the 30°-Stokes' wedge centered about -60° (see Fig.[2|. 
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FIG. 1: Stokes' wedges (shaded regions) in the complex-ip 
plane in which the integration contour C for the integral 
J c dip exp (— <y3 2 ) terminates. This integral does not exist if 
C terminates in an unshaded wedge. 




FIG. 2: Stokes' wedges (shaded regions) of angular open- 
ing 30° in which the integration contour C for the integral 
J c dp exp (—p 6 ) may terminate. The integral has two possi- 
ble real values, one for which the contour connects the pair of 
wedges centered about the real axis and the other for which 
the contour connects the lower pair of wedges. 

The contour C for J c dtp exp (— <p 6 ) must join a pair 
of 'PT-symmetric Stokes' wedges (wedges that are sym- 
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FIG. 3: Stokes' wedges (shaded regions) in which the integra- 
tion contour C for J c dtp exp (— <£> 4n+2 ) may terminate. When 
n — 2 (left figure) there are three pairs of Stokes' wedges and 
when n = 4 there are four pairs of Stokes' wedges. 



metric about the imaginary axis) or else the integral is 
not real. A third pair of 30°-Stokes' wedges, one cen- 
tered about 60° and the other centered about 120°, are 
not shown in Fig. [5] the integral exists if the contour G 
connects this pair of Stokes' wedges, but this case is not 
new; it is just the complex conjugate of the configuration 
in which C connects the —120° and —60° wedges. 

The cases n = 2 (three pairs of 18° Stokes' wedges) 
and n = 3 (four pairs of 12.8° Stokes' wedges) are shown 
in Fig. [3] In the former case the J c dip exp (— ip 4n+2 ^j has 
three independent real values; in the latter case it has 
four independent real values. 

Returning to the quantum field theory with vacuum 
persistence function given in ([3]), we vary the action in 
the exponent and obtain the Euclidean field equation in 
the presence of the external c-number source J(x): 



V 2 <Kz) + g[tf>(x 



l4n+l 



J(x) 



(4) 



This field equation is unique; it does not depend on the 
choice of complex contour G. 

The expectation value of Q in the vacuum state is 

- V 2 Gr(:r) + g([<t>(x)] in+1 ) /Z[J] = J(x), (5) 

where G\ [x) is the connected one-point Green's function: 

8lnZ[J] (<t>(x)) 



SJ(x) 



Z[J] 



T>4> <p{x)e 



(6) 



This expectation value depends on the choice of metric, 
but in Ref. [10] it is shown that the path integral auto- 
matically gives the expectation value with the appropri- 
ate metric. Thus, if the integration contour terminates 
in the wedges containing the real axis, then the metric 



uses the conventional Dirac adjoint f, and if the contour 
terminates in another pair of Stokes' wedges, then the 
metric uses the CVT adjoint of the corresponding non- 
Dirac-Hermitian "PT-symmetric field theory 

To derive the Dyson-Schwinger equations for the con- 
nected Green's functions of the quantum field theory, we 
express the second term on the left side of ^ in terms of 
the higher connected Green's functions. The technique is 
standard (see, for example, Ref. [12]); one differentiates 
repeatedly with respect to the external source J(x) and 
uses the formula for the rt-point Green's function in the 
presence of the external source J: 

G n {x, y,z,...) = S n /[SJ(x)SJ(y)SJ(z) • ■ ■ ] In Z[J}. (7) 

We must truncate the Dyson-Schwinger equations in 
order to obtain a closed system. We consider here just 
the first two equations and neglect contributions from 
Green's functions beyond G^^, y)- This truncation gives 
the mean-field (or one-pole) approximation to the two- 
point Green's function. (Including higher Green's func- 
tions does not change any qualitative conclusions of this 
paper.) Thus, we repeatedly differentiate with respect to 
J(x) and get the following sequence of equations: 

(1) = Z[J], {cj>{x))=G l {x)Z[Jl 
([4>(x)} 2 ) = ([G 1 (x)] 2 + G 2 (x,x))Z[J], 
([<P(x)} 3 ) = ([Gr(x)] 3 +3G 1 (x)G 2 (x,x)) Z[J], 
MX)] 4 ) = ([Gr^ + eldOE)] 2 ^,*) 

+-i[G 2 (x,x)] 2 )Z[J], 
([cj>(x)} 5 ) = ([Gr(x)] 5 + 10[G 1 (x)] 3 G 2 (x,x) 

+15G 1 [G 2 (x,x)} 2 )Z[J}. (8) 
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These expressions have a simple form as polynomials 
P n (t) in the variable t = G\{x) J ^/G^ix, x), 



{[ct>{x)] n )^[G 2 (x,x)] n ' 2 Z[J]P n {t), 



(9) 



where P n (t) = (— i) n He n (it) are Hermite polynomials of 
imaginary argument: Po(t) = 1, P\(t) = t, P 2 (t) = t 2 + l, 
P 3 (t) = t 3 +3t, Pi{t) = < 4 +6< 2 +3, P 5 (t) = t 5 + 10t 3 + 15i. 
Next, we insert ^ into ^ and obtain 

- V 2 Gi(x) -^[G 2 (x,:r)] 2n+1 / 2 Hc4„ + i^) = J(x). (10) 

At J = translation invariance is restored, and G\(x) 
and G 2 (x, x) become the numbers G\ and G 2 (0). Thus, 
the first of the truncated Dyson-Schwinger equations is 



He 



4n+l 



0. 



(11) 



Note that the argument of He4„+i remains invariant if 
wave-function renormalization is performed. 

To obtain the second Dyson-Schwinger equation we dif- 
ferentiate (llO]) with respect to J(y) and set J = 0: 



(-V 2 + M 2 )G 2 (x-y) = 6 D (x-y), (12) 
where the renormalized mass is given by 

1 (13) 



M 2 = [G 2 (0)] 2n He^ +1 iG 1 /^G 2 {Q) 



We solve ( 11 )— (|13|) simultaneously: First, we Fourier 



transform |12j) and find that in /^-dimensional Euclidean 
space G 2 (p) = l/(p 2 + M 2 ). Thus, for < D < 2 we get 
the finite result G 2 (0) = M D - 2 r(l - D /2)2- d tt- d / 2 . 
Second, we note that the Hermite polynomial He4„ + i is 
odd and only has real roots. There are two cases: Either 
(i) Gi = 0, which is the conventional Dirac-Hermitian 
parity-invariant solution to the Dyson-Schwinger equa- 
tions, or (ii) we get 4n new parity-violating nonzero val- 
ues for the one-point Green's function: 



(14) 



where the dimensionless number rj (J = 1, . . . 2n) is one 
of the 2n positive roots of He4„+i. Finally, we use the 



identity He 4ra+1 = (4n + l)He4„ in (13) and use the 



interlacing-of-zeros property of the Hermite polynomials 
to prove that there are exactly n new positive values of 
M 2 corresponding to the nonzero values of G\j. This 
demonstrates the connection between pairs of Stokes' 
wedges and solutions to the Dyson-Schwinger equations. 

For example, when D — 1 in a 4> 6 model, r$ = and 
n = 2.85697, and corresponding to these roots the di- 
mensionless renormalized masses are M — 1.39158 and 
M = 2.25399. Thus, there are two families of particles: 
One particle (associated with a nonvanishing G^i) has 
a mass 1.62 times larger than that of the other particle 



(associated with a vanishing Gi,o)- This ratio increases 
rapidly as a function of the space-time dimension D; for 
example, for D = 0.0, 0.5, 1.0, 1.5, 2.0, 2.5 this ratio 
takes the values 1.38, 1.47, 1.62, 1.90, 2.62, 6.88. 

To conclude, while a flavor symmetry group is conven- 
tionally introduced to describe families of particles, we 
have shown that such families can arise naturally from 
the monodromy structure in the complex-field plane as- 
sociated with rotation from one Stokes' wedge to another. 
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